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The dynamics of non-spherical rigid particles immersed in an axisymmetric random flow 
is studied analytically. The motion of the particles is described by Jeffery's equation; the 
random flow is Gaussian and has short correlation time. The stationary probability den- 
sity function of orientations is calculated exactly. Four regimes are identified depending 
on the statistical anisotropy of the flow and on the geometrical shape of the particle. 
If A is the axis of symmetry of the flow, the four regimes are: rotation about A, tum- 
bling motion between A and —A, combination of rotation and tumbling, and preferential 
alignment with a direction oblique to A. 



1. Introduction 

Non-spherical solid particles suspended in a moving fluid rotate and orient themselves 
under the action of the velocity gradient. Even at low concentrations, the orientational dy- 
namics of non-spherical particles can influence the rheological pr operties of a suspension, 
name ly the intrinsic viscosity and the normal stress coefficients ( Bird et al. 1977t L arson! 
1999). This phenomenon has diverse practical applications. In the turbulent regime, for 
instance, the injection of rodlike polymers in a Newtonian fluid can produce a consider- 
able reduction of the turb ulent drag, with this effec t being routinely exploited to reduce 
energy losses in pipelines ( Gvr fc Bewersdorfj 1995 ). The study of the orientation of par- 
ticles immersed in a fluid also has numerous applications in the natural sciences. Amongst 
them it is worth mentioning the swimming motion of c ertain biological micro-organisms 
(jSaintillan fc Shelley 2008: Koch fc Subramanianll201ll ) and the formation of ice crystals 
in clouds ( Chen fc Lamb|[l994) . This latter phenomenon plays a crucial role in processes 
such as rain initiation and radiative transfer. 

The starting point for understanding the properties of a dilute suspension is the mo- 
tion of an isolated particle in a given flow field. Analytical results on the dynamics of 
a non-sph erical particle have been obtained for various laminar flows, both steady and 
unsteady. Ijefferv ( 1922 ) derived the equations of motion for an inertialess ellipsoid in a 
steady uniform shear flow at low Reynolds number. For a spheroid (i.e. an ellipsoid of 
revolution), Jeffery showed th at the axis of sym metry of the particle performs a peri- 
odic motion on a closed orbit. Brethertonl (|l962h subsequently extended Jeffery's anal- 
ysis to particles of a more general shape: he demonstrated that, except for certain very 
long particles, the dynamics of any body of revolution transported by a low-Reynolds- 
number shear is equivalent to that of a spheroid with an effective aspect ratio. Amongst 
bodies of revolution, rigid dumbbells and rods hav e received a systematic investigation 
(|Bird. Warner fc Evanslll97ll: IDoi fc Edwardslll986T) . 

In Jeffery's il922j) and Bretherton's (| 1962T) derivations, not only fluid and particle 
inertia are disregarded, but also Brownian fluctuations due to the collisions of the particle 
with the molecules of the fluid. Nevertheless, if a particle is sufficie ntly small, mo l ecular 
diffusion does influence its orientational dynamics, as was shown bv lLeal fc Hmchl (|l97lh 
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and Hinch fe Leall (1972). The review article bv I Brenner (1974) collects analytical results 
on the motion of rigid neutrally buoyant bodies of revolution subject to a uniform velocity 
gradient and to Brownian fluctuations. A more recent review on this problem can be 
found in IPetrid (|l999h . Finally, even in simple laminar flow s, the orien tation of particles 
can form complex patterns; this behaviour was studied bv lSzeril(|l993h in the context of 
the theory of dynamical systems. 

In turbulent flows and in some chaotic flows, the velocity gradient exhibits small-scale 
fluctuations. Thus, depending on the size of the particles, Jeffery's assumption of a spa- 
tially uniform velocity gradient may not be applicable. To measure the probability of 
orientations, experiments have generally used relatively large particl es, and therefore the 
velocity gradient could vary appreciably over the size of a particle (e.g. Krushkal fe Gallfly! 
19881: iBernstein fe Shapirolll994l:lNewsom fc Bry Icll998tlPa rsheh. Brown fe Aidunl l2'o05 



Parsa et al. 2011 ; Zimmermann et al. 1201 lal bh . Accordingly, theoretical studies were 



mainly concerned with the derivation of model equations for the probability of ori- 
entations, in which turbulent fluctuations were treat ed as an effective isotropic diffu- 
sion term ( Olson fe Kerekes 1998 ; Shin fc Koch 20051 ). Jeffery's approach, however, re- 
mains applicable to chaotic or turbulent flows provided that the particles are sufficiently 
small. The orient ation dynamics of tiny elongated particle s was studied numerically 
in channel flows dZhang et a/.ll200lt iMortensen et qLI l2008allbh . in isotropic turbulence 



(IShin fc Kochll2 005; Pumir & Wilkin sonll2011tlParsa et aLll2012D . and in chaotic velocity 
fields (|Wilkinson. Bezuglvv fc MehlisJl2009l) . In contrast with the case of laminar flows, 
few analytical results seem to exist for the probability distribution of orientations of 
small particles transported by a turbulent or chaotic flow. iTuritsvnl (|2007l ) examined the 
tumbling motion of rodlike polymers in a random flow res ulting from the superposition 
of a mean shear and of white- in-time isotropic fluctuations. I Wilkinson fc Kennardl (|2012l) 
recently studied the alignment of rods with vorticity in a turbulent isotropic flow. 

Here, the probability density function (p.d.f.) of orientations is derived exactly for small 
particles transported by a random flow with axisymmetric statistics. The particles are 
general bodies of revolution possessing fore-aft symmetry. The axisymmetry of the flow 
means that the velocity field is statistically invariant under arbitrary rotations of the axes 
about a given di rection as well a s under reflections in p lanes containing that direction 
or normal to it dBatchelorl 1 1946 \_ Chandrasekharljl950[) . Axisymmetry is the simplest 
form of statistical anisotropy (jBiferale fc Procaccial 2005L Chang. Bewley fc Bodenschatz 
2012 ). and is found in rotating, stratified, or wind-tunnel turbulence ( Lindborgl Il995l ) . 
Furthermore, the random flow is assumed to be Gaussian and to have zero correlation 
time. The assumption of temporal decorrelation is adequate when the correlation time 
of the flow is short compared to the characteristic time scale of material-line-element 
stretching. This assumption, albeit restrictive, allows a fully analytical solution of the 
problem. 

The evolution equation for the orientation vector of a non-spherical rigid particle is 
introduced in § [2] Section [3] is devoted to the derivation of the Fokker-Planck equation 
for the p.d.f. of the orientation angle. Its stationary solution is studied in § 2] Some 
conclusions are drawn in § [5] 



2. Orientation dynamics 

The particles considered here are rigid bodies of revolution possessing fore-aft symme- 
try (although in the literature such particles are commonly referred to as "axisymmetric" , 
this terminology will be avoided here not to generate confusion; the term "axisymme- 
try" will be reserved to the statistical invariance of the velocity field). The particles 
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are of uniform composition and are suspended in a Newtonian fluid of the same den- 
sity. Furthermore, the inertia of the particles as well as hydrodynamic particle-particle 
interactions are disregarded, and no externally imposed force or couple influences the 
dynamics. In particular, it is appropriate to disregard hydrodynamic interactions when 
the suspension is sufficiently dilute. 

The undisturbed motion of the fluid is described by the velocity field v(x,t). The 
size of the particles is assumed to be small compared to the typical length over which 
the velocity gradient Vi? = (3jUj)i<i )i j<3 changes (3^ = d/dxj). For turbulent flows, 
this assumption is satisfied if the particles are smaller than the Kolmogorov scale, where 
viscosity smooths out the velocity fluctuations; the Reynolds number of the flow at the 
scale of a particle is thus less than 1. Given their small size, the particles also experience 
Brownian collisions with the molecules of the fluid. 

In a sufficiently dilute suspension, attention can be restricted to the dynamics of a 
single isolated particle. The configuration of a body of revolution is determined by the 
position of its centre of mass, r c (t), and by the orientation of its axis of revolution, which 
is specified by a unit vector N(t) parallel to the axis itself. As the particles are subject to 
Brownian fluctuations, their dynamics is random even in a laminar flow. Consider first 
a deterministic velocity field or a given realisation of a random velocity field. On the 
above assumptions, t he centre of mass moves according to the following equation (e.g. 
iDoi k Edwardlll986l ): 

r c (t) =v(r c (t),t) + v^rCW. I 2 - 1 ) 
where T>t > is the translational diffusion coefficient and is three-dimensional white 
noise, i.e. a Gaussian stochastic process with 

(C(*))=0 and (Q(t + T)( j (t)) = 6 ij 5(T) (2.2) 

for all t,r > and i,j — 1,2,3. The orientation vector satisfies the following stochastic 
differential equation (summation over repeated indexes is implied): 

N> = KijWNj - K pq (t)j^-Ni + y/v^SaiN) o &(t), |JV(0)| = 1, (2.3) 

where 

K(t)=Q{t)+jE{t) (2.4) 

with 

m - M^w, E(t) . + (2.5, 

and G(t) = V«(r c (t),t). Thus, Q(t) and E(t) are the vorticity tensor and the rate-of- 
strain tensor evaluated at r c (i) and n(t) is an effective Lagrangian velocity gradient. 
The scalar constant 7 depends on the geometrical shape of the particle. For |-y | < 1, 
the evolution equatio n can be mapped into that of a spheroid with aspect ratio equal 



to + 7)7(1 — 7) (jBrethertonl Il962f ) . Prolate spheroids are obtained for < 7 < 1, 
oblate spheroids for —1 < 7 < 0. Special cases are: spheres (7 = 0), rigi d dumbbells 
(7 = 1), rods (7 = 1), and disks (7 = —1). Furthermore, it was shown bv IBretherton 



(1962) ;hat in principle there exist very long particles for which \j\ > 1. In (|2.3p . the 



random vector £(t) is three-dimensional white noise and hence has the same properties 
as C(£), but is statistically independent of it. The matrix Z(n) has the following form: 

r(n) = / -nn/\n\ 2 (2.6) 



and Vr > is the rotary diffusion coefficient. By using the Cauchy-Schwarz inequality, 
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it is easy to check that E(n) is positive semi-definite. The symbol o indicates that the 
stochastic term in ()2.3[) is understood in the Stratonovich sense. 

Equation (|2.3I) is Jeffery's equation for the orientation vector of a body of revolution 
with the addition of a stochastic term modelling Brownian fluctuations. The stochas- 
tic term is chosen in such a way as to produce isotropic diffusion of N(t) on the unit 
sphere, so that |iV(f)| is preserved in time (see appendix IA1 for more details). It is worth 
remarking that a Brownian term of th e same form has been used to model the turbulent 
fluctuations of the velocity gradient ( Krushkal fc Gallfly 1988t Olson fe Kereked 1998 ; 
Shin fc Kochll2005l) or to descri be particle-particle interactions b oth in semi-dilute sus- 



pensions (Doi fc Edwardslll986l ) and in concentrated suspensions ( Doi fc Edwards! 1978 ; 
Kuzuu fc Doilll980h . 



Equation (12. 3|) can be generalised to the case of a homogeneous axisymmetric random 
flow. The velocity field transporting the particle is Gaussian and has zero mean and 
correlation: 

{v i {x + r,t + r)v j (x,t))=Q ij (r)5(T), i,j = 1,2,3. (2.7) 
The form of the correlation guarantees that v(x,t) is statistically homogeneous in space. 
Additionally, the velocity field is assumed to be incompressible (V • v = 0) and statis- 
tically axisymmetric with respect to the direction specified by the unit vector A. The 
tensor Q(r) must then take the form: 

Qij{r) = Anrj + BSij + CX.Xj + D(nXj + A^), (2.8) 

where A, B, C, D are smooth fu nctions of |r| 2 and (r • A); A, B, a nd C are even 
in (r • A), while D is odd in (r • A) ( Batchelor 1946t Chandrasekhar 1950h . Furthermore, 
the functions A, B, C, D are not inde pendent and satisfy certain differential relations 
( Batchelorlll94d : Chandrasekharl Il950l ) . The velocity field defi ned above i s an a xisym- 
metric generalisation of the isotropic random flow introduced bv lKraichnar] (|l968l ) in the 
co ntext of passive turbule nt transport. The same axisymmetric velocity field was used 
bv lShaqfeh fc Ko ch (1992) to study polymer stretching in flows through random beds of 
fibres. 

The velocity gradient is also Gaussian and zero-mean; the single-point two-time cor- 
relation can be derived by using the statistical homogeneity of the velocity field: 

(djVi(x,t + T)dgV p (x,t)) = r ijpq S(T) i,j,p,q =1,2,3 



with 



A 



3 Q v 



drjdr q 



(2.9) 



(2.10) 



(2.11) 



Substituting ff2T8]> in (|2~10| yields (see equation (5.12) in iBatcheloil (|l946h ): 

Aj>9 =(d + ^ajdjqSip - a(S pq Sij + S iq 5 P j) + (b + c + 5d)5 jq XiX p - b5 ip XjX q 
— d[(5 pq Xi + 5i q X p )Xj + (S P jXi + SijX p )X q ] — cXiX p Xj X q , 

where a, 6, c, d are real constants. For the sake of simplicity, A is taken in the direction 
of the third axis, i.e. A = (0,0, 1). Then, the coefficients in Bij pq arc written: 

2a = A212 _ Am d = ^3333 — Am 

(2.12) 

b = A212 - A313 2a + 6 + c + 4d = A131 ~ A333- 

The two-time correlation of the co mponents of the vorticity u> = V x v can be expressed 
in terms of a, c, and d as follows (|Batchelorlll946l p. 490): 

(u)i(x,t + T)u3i(x,t)) = (u 2 (x,t + T)u2(x,t)) = (10a + c + 9d)S(r) (2.13) 
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and 

(w 3 (x, t + t)oj 3 (x, t)) = (10a + 2d)*(r). (2.14) 

The coefficients a, 6, c, d are not free; they are constrained by the following inequalities 
(see appendix IB]) : 

ja + d>0 ifa>0, 
[5a + d>0 ifa<0, 

4a-6 + d>0, 4a + o + c + 6d>0, (2.156) 
15a 2 - 6 2 - (6 - d)(c + 5d) + 2a(2c + 13d) > 0. (2.15c) 

If a > and b = c = d = 0, t hen Fj^g gives t he single-point correlation of the gradient 
of an isotropic velocity field ([Robertson 1940). Thus, a determines the intensity of the 



isotropic part of the gradient and 6, c, d control the statistical anisotropy of the flow. 

The orientation dynamics of the particle depends on the velocity gradient evaluated 
at r c (t), which was denoted as G(t) in § [2] In virtue of the (^-correlation in time and the 
statistical homog eneity of the flow, G(t) has the same te mporal statistics as Vw(a5,t) 
for any given x (|Falkovich. Gawedzki fc Vergassolall200ll ) (the presence of white noise 



in (|2.1[) does not modify the statistics of G(t)). The components of K(t), defined in (I2.4[) . 
are a linear combination of the components of G(t), and consequently K,(t) is a Gaussian 
process with zero mean and correlation: 

{nij(t + T)n pq (t)) = K lJpq 8{T) (2-16) 

with 

Kijpq — ^ [-^ijpg ^ij<lP ^jipq~^~-^jiqp~^~'^'~f{-^ijpq ^jiqp) 1 {J^ijpc[~^~ ^ijqp~^~ ^jipq~\~ ■^jiqp)\ t 

(2.17) 

The form of Kij pq can be derived by substituting (|2.4p and (|2.5[) into the left-hand 
side of (|2.16[) and by using (|2.9p . The tensor Kij Pq is positive semi-definite, for it is 
the covariance of a second-order Gaussian tensor. Thus, in the evolution equation for the 
orientation vector, n(t) plays the role of a multiplicative tensorial white noise. As n(i) can 
be thought of as an approximation of a real noise process in the limit of zero correlation 
time, the corresponding terms in (|2.3[) must be interpreted in the Stratonovich sense (e.g. 



Kloeden fc Platenlll992l p. 227). The stochastic differential equation for the orientation 



vector of a non-spherical particle can then be rewritten as follows: 

N t = M ipq {N) o Kpq (t) + y/v^Ei^N) o &(t), (2.18) 

where Mi pq (n) = (Si P Sjk — Sik8jp)njnkn q /\n\ 2 and the initial condition iV(0) is such 
that |JV(0)j = 1. 



3. Fokker Planck equation for the probability density function of the 
orientation angle 

As v(x,t) is statistically invariant under spatial translations, the p.d.f. of N(t) taking 
the value n = (7^1,712,713) at time t is independent of r c and is thus denoted by f(n;t). 
The Ito equation equivalent to (I2.18[) is 

Ni = ft(JV) + M ipq {N)K pq (t) + s/V^EaiN)^) (3.1) 
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Figure 1. Orientation of a non-spherical particle. 



with 



1 3 

Pd n ) = 2 K mn P qMj pq (n)—M imn (n) + 



^ Sjk{ n)±-S ik {n). 



Consequently, f(n; t) satisfies the Fokker-Planck equation: 



dt 



orii 2 oniOrij 



>«(")/] 



with 



a i:j {n) 



in 

{n)M Jpq (n) +T> R E ik (n)E jk (n) 

ill 

(n)Mj P q(n) +U R 2J ij (n), 



(3.2) 



(3.3) 



(3.4) 



where the last equality follows from (|A 31) . Equations fl3.1[) and Q3.3P can be derived 
from p,18p by using the formal rules Kij(t)dt = 0(^fdi) and Kg (t)dt K va (t)dt = Kij pq dt 



and by proc eeding as in the case of a vectorial white noise (see lGardinerlll983l and the 
appendix in lFalkovich et cd. 2001 ). The diffusion tensor a. is positive semi-definite as a 
consequence of the positive semi-definiteness of AC and E. 

To study the orientation dynamics of a non-spherical particle, it is convenient to move 
from Cartesian coordinates (n\, ri2, n^) to spherical coordinates {n, i9, tp) according to the 
usual transformations: 



'3- 



d = arctan 



(\/ n i n i/ n 3) ! f — arctan(ri2/wi) (3.5) 



with ^ n, ^ d ^ tt, and ^ <p < 2ir (figure [T]). On account of the fixed length of the 
orientation vector, the probability density function of orientations must take the form 
f(n, i?, tp; t) = t)8(n — L) with L = 1. Thus, ip^i f\ t) sm$ dz? dp is the probability 

of the particle being oriented at time t within an elementary solid angle sin t9 dD dp 
of ($,</?). In addition, the following normalisation holds: 



'2 7T 



ip{"d, ip; t) sin $ d-d dp = 1 



Vi > 0. 



(3.6) 



The function W{d : p\t) = ip("d, <p; t) sini? satisfies a new Fokker-Planck equation, which 
can be derived from (|3. 3|) by using the tr ansformation formulae for the drift and diffusion 
coefficients under a change of variables ( Riskenl 19891 p. 88), by writing f(n,'&,(p;t) = 
tp(i!},<p;t)5(n — L) with L = 1, and by integrating the resulting equation with respect 
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to n. The final result is: 



f = + + (3-7) 



where 



B*(fi) = Pi(n) £- + = [ Ml + M2 sin 2 (tf) + M3 sin 4 (d)] cot tf, (3.8a) 

9$ dd 

Aw>{-&) = a tJ (n) — — = 2/x a + ^4 sin 2 (tf) + /i 3 sin 4 (tf), (3.86) 
on, cm^ 



.4^,(0) = ay (n) |^ J^- = /is + 2/xi csc 2 (tf) (3.8c) 
an, onj 



with 



/H = i [2(5 + 3 7 2 )a - 4 7 6 + (1 - 7 ) 2 c + (9 - IO7 + 5 7 2 )d + 4T>' R ] , (3.9a) 

,j 2 = f 7 [26+(l- 7 )c+(5- 7 )d], (3.96) 

M3=c 7 2 , (3.9c) 

(H = 7 [26+ (1 - 7 )c+5d], (3.9d) 

M5 = \ [46 7 - (1 - 7 ) 2 c - (1 - 7 )(7 - 3 7 )d] , (3.9e) 

and 2?^. = T>r/L 2 with L = 1 (the numerical values of T>' R and Pjj coincide, but their 
physical dimensions are different). Equations (|3.8[) involve the Jacobian and the Hessian 
of the transformation from Cartesian to spherical coordinates, which can be calculated 
from (|3.5j) . Inequalities (|2.15j) guarantee that Am^) and A vv {{f) are strictly positive 
(appendix |B|). 

The contribution to (|3.7j) due to the isotropic part of f is a Laplace-Beltrami term 
with diffusion coefficient proportional to a; this contribution is of the same form as that 
coming from £(t). The isotropic component of the flow and the Brownian fluctuations 
therefore have the same effect on the orientation statistics of the particle. 

Since 1? and ip are angular variables, the boundary conditions for ^(d, ip; t) are periodic: 

<P(i?, <p; t) = + 2n, cp; t) and &(<&, ip; t) = <p + 2n; t) (3.10) 

for all i9, ip, and t. The long-time properties of </j;i) can be deduced from (|3.7j) . If 
the partial derivative with respect to time is dropped, then (|3.7j) is invariant under the 
transformations d «-» 27t — $ (reflections with respect to planes containing A) and ■& 4-» 
7t — $ (reflections with respect to planes orthogonal to A). Moreover, the coefficients B$, 
Atftf , and A Vip do not depend on ip (invariance under rotations about A) . These properties 
of (|3.7|) are a natural consequence of the statistical axisymmetry of the velocity field and 
translate into analogous properties of the stationary p.d.f. of orientations. 

The invariance of (|3.7p under rotations about A can be used to derive a one-dimensional 
Fokker-Planck equation for the marginal p.d.f.: ^(■d;t) = ?/>(#;£) sin # with il)(d]t) = 
J ipifl, (p; t)d<p. Integrating (|3.7I) with respect to <p from to 2n and making use of (|3. 10[) 
yield: 

d$ _ 3 

~a7 ~~ 



1 9 2 



Mm +29^ AmW . (3.11) 



The solution of the above equation must be normalised and periodic: t) = 'P(i3+2n; t) 
for all #, t. A direct consequence of (|3.11j) is that, along the trajectory of the particle, the 
time evolution of #(£) is decoupled from that of <p(t) and is described by the stochastic 
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ordinary differential equation: 



0(4) - B*(0(t)) + VA?tf(0(4))?7(4), 
where 77(4) is white noise. 



< 0(4) < 7t, 



(3.12) 



4. Stationary statistics of the orientation angle 

It was argued in §[3] that the stationary solution of (|3.7[) . ^ st , does not depend on tp. 
Therefore, solves the equation: 

\ ^ [AM^st] - ^ [B*(W«] = 0- (4-1) 

Thanks to the periodic boundary conditions p.lOp . !Pst(0) takes the following simple 
form (see appendix ICj): 



A/' 



>W(0) 



exp 



gg(g) 
# Am{z) 



where 0o ^ 7t and A/" is a normalisation constant such that 

27T / ^ st (0)d0 = 1- 

Jo 



(4.2) 



(4.3) 



In (|4.2p . the choice of 0o is in fact unimportant, since it only modifies the normalisation 
constant. 

Only the case /13 7^ is considered here; the case /13 = is examined in appendix [Pi 
even though no new physical regimes emerge in this latter case. For ^3 ^ 0, the integral 
in (|4.2I) can be calculated by using the change of variab le y = sin 2 (a;) and formulae 2.172, 
2.175(1), and 2.177(1) of iGradshtevn fc Rvzhikl (jl965l ). The final result is: 



0st(0) = 



where 



exp 



fj,2 A*4 arctan 



A^(0) 

/i 4 + 2fi 3 sin 2 (0) 



(4.4) 



A 



xW = { exp 



^4 



M2 



/i4 + 2(i 3 sin (0) 



/ =: 2\ + ^4 + 2^3 sin 2 (0) 



-A - ^ 4 - 2^3 shr(0) 



(z4>0) 



(z4 = 0) 



(z4<0) 



(4.5) 



with A = 8/ii/i3 — It is shown in appendix [Bl that "081(0) is bounded for all values 
of a, b, c, d, and 7. The stationary p.d.f. of orientations satisfies tfj s t(2n — 0) = ip s t($) 
and i/j a t(n— 0) = V>st(0)- These properties are a consequence of the statistical symmetries 
of the carrier flow, as was noted after (|3.10[) . 

Since A can be written as A — a 2 F(-j,b/a,c/a,d/a,T>' R /a), Vst(0) only depends on 
the ratios b/a, c/a, d/a, T>' R j a (and on 7). The same conclusion could have been reached 
by rescaling 4 by a -1 in (|3.7I) . 

For a spherical particle (i.e. 7 = 0), ^ 2 and ^4 vanish. Consequently, x(0) = 1, 
-4i?i?(0) = const., and hence V>st(0) = (47t) _1 , in accordance with the fact that all orien- 
tations are equally probable for a sphere. Similarly, if a or V R are much greater than 6, 
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Figure 2. Stationary p.d.f. of the orientation angle (a) Rotation about A: b/a = 4.8, c/a = 55, 
d/a = 0.85, V R /a = 10~ 2 , and 7 = -1 (solid line), 7 = -0.75 (dashed line), 7 = -0.5 
(dot-dashed line), 7 = —0.25 (dotted line), (b) Tumbling motion between A and —A: b/a — 4.75, 
c/a = 10, d/a = 0.9, T>' R /a = 10~ 2 , and 7 = 1 (solid line), 7 = 0.75 (dashed line), 7 = 0.5 
(dot-dashed line), 7 = 0.25 (dotted line). The normalisation coefficient J\f has been computed 
numerically according to (|4.3|l . 



c, d, i.e. if the isotropic component of the velocity field or the Brownian fluctuations 
prevail on the anisotropic component of the flow, then fi% ^> fa, i = 2, . . . , 5 and V'st(^) 
weakly depends on d regardless of the shape of the particle. In the following, therefore, 
7 is assumed to be nonzero and b, c, d are of the same order of magnitude as a and T>' R 
or greater. 

The behaviour of ip s t($) can be deduced from that of its first derivative. By us- 
ing Ami'ff) > 0, it can be shown that, for all values of A, 

-j^stW = 7sin(2i?)[c 7 cos(2i?) - tr]h{€), (4.6) 

where 

cr = 26 + c+(5 + 7)d (4.7) 

and is a strictly positive function for all ^ d ^ 7t. Four regimes are then identified 
depending on the properties of the axisymmetric flow and on the geometrical shape of 
the particle: 

(i) Rotation about the axis of symmetry of the flow. For the following values of the 
parameters: 

|er| > |c7| and 70- < 0, (4.8) 

the function tp st only has three extrema in $ = 0, 7t/2, 7t (indeed the equation c 7 cos(2i9) = 
a has no solution). More precisely, "0st(^) has a maximum in 7t/2 and two minima in 
and 7t (figure^. Thus, in this regime the particle rotates about the direction A; the level 
of alignment with the plane perpendicular to A decreases as the degree of anisotropy of 
the flow vanishes or the shape of the particle approaches the spherical one (figure [2J . 

(ii) Tumbling motion. In the following regime: 

\a\ > \cj\ and 70- > 0, (4.9) 



ipst has three extrema: a minimum in 7t/2 and two maxima in and in n (figure [2]). 




°0 nil n °0 7i/2 n 

e e 

Figure 3. Stationary p.d.f. of the orientation angle (a) Preferential alignment with a direction 
oblique to A: 7 = 1, T>' R /a — 10 -2 , b/a = —6, d/a — —0.9, and c/a = 55 (dashed line), c/a — 25 
(solid line), c/a — 11 (dot-dashed line), (b) Combination of rotation and tumbling: 7 = 1, 
T>' R /a = 10~ , c/a = —10, d/a = 1, and b/a — 1.1 (dot-dashed line), 6/a = 2 (solid line), 
fo/a = 3 (dashed line). The p.d.f. has been normalised numerically according to (14. 3|) . 



The particle tumbles between the direction parallel to A and that antiparallel to A. The 
probability of the orientation angle $ being in the neighbourhood of or tc depends on 
the anisotropy degree of the flow and on the shape of the particle (figure [2J . 

(iii) Preferential alignment with a direction oblique to the axis of symmetry of the flow. 

If 

|er| < \cr/\ and c> 0, (4.10) 

then ip si has three minima in 0,7t/2,7t and two maxima in and re — where < 
t?* < 7t/2 is such that 




sin^ = W- 1-- . (4.11) 



The particle therefore spends most of the time at an angle (or 71 — with respect 
to A (figure [3]). 

(iv) Combination of rotation and tumbling. For 

H < |c 7 | and c < 0, (4.12) 

the function tp s t has three maxima in 0, 7t/2, tc and two minima in and re — 19* with 
denned in (|4.11j) . In this regime, the particle preferentially lies either in the plane per- 
pendicular to A, in the direction parallel to A, or in the direction antiparallel to it. 

Naturally, the coefficient a controlling the intensity of the isotropic component of the 
flow does not play any role in the above classification; the dynamical regime is selected 
by b, c, d, and by the shape coefficient 7. Whereas sufficiently elongated or flattened 
spheroids can be strongly aligned in regimes (i) and (ii), the ability of the flow to orient 
particles is weaker in regimes (iii) and (iv). I n these regimes, st ronger alignment can be 
obtained for |7| > 1 (figure 0]). Nevertheless. iBrethertonl (1962) observed that particles 



with I7I > 1 may be unrealistic, albeit conceivable from a purely geometrical point of 
view. 

The form of V'st(^) simplifies considerably if (divi(x,t + r)divi(x,t)) — (33^3(0;, t + 
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Figure 4. Stationary p.d.f. of the orientation angle $ for |^| > 1. (a) Preferential alignment 
with a direction oblique to A: 7 = 5, T>' R /a = 10~ 2 , b/a = —6, c/a — 25, d/a = —0.9. (6) 
Combination of rotation and tumbling: 7 = 5, T>' R /a = 10 -2 , b/a = 2, c/a = —10, d/a — 1. 



t)9 3 v 3 (:c, £)), i.e. if d = 0. In this case, ^ 2 = 3^/4 and hence VstW = •A/"[.4#,9(i?)] _3//4 ; 
furthermore, cr does not depend on 7. It then follows from (|4.8|l and (14.91) that if the 
particles with shape coefficient 7 rotate in the plane orthogonal to A (resp. tumble), the 
particles with shape coefficient —7 tumble (resp. rotate in the plane orthogonal to A). 
By contrast, regimes (iii) and (iv) are independent on the sign of 7, i.e. on whether 
the particle is elongated or flattened, although changes when the sign of 7 changes. 
Moreover, if the particles with shape coefficient 7 rotate about A (resp. tumbles), then 
all particles with \j\ < \j\ and sgn(7) = sgn(7) rotate about A (resp. tumble). Similarly, 
if the particles with shape coefficient 7 are in regime (iii) (resp. in regime (iv)), then 
all particles with I7I > I7I are in regime (iii) (resp. in regime (iv)). These properties, 
however, do not generally hold true if d is nonzero. For instance, if b/a = —4.4, c/a = 0.5, 
d/a = 1.5, then particles tumble for 1 J? 7 > 0.8, they have a preferential orientation 
for 0.8 > 7 > 0.4, they rotate for 0.4 > 7 > 0, and they tumble again for > 7 ^ — 1. 



5. Conclusions 

Axisymmetric turbulence arises as one of the simplest frameworks in which to study the 
orientation dynamics of non-spherical particles. On the assumptions of Gaussianity and 
short correlation in time, it was shown analytically that the dynamics of a non-spherical 
particle immersed in a random axisymmetric flow exhibits four regimes: rotation around 
the axis of symmetry of the flow, tumbling, combination of rotation and tumbling, and 
preferential alignment with a direction oblique to the axis of symmetry of the flow. The 
regime is selected by the form of the anisotropic component of the flow and by the 
geometrical shape of the particle. If the flow is weakly anisotropic or if the particle is 
almost spherical, the mathematical description of ip s t (i?) in terms of minima and maxima 
remains formally valid, but the above physical classification loses its meaning, since ipst ($) 
does not differ appreciably from the uniform distribution. 

The tumbling motion of a non-spherical particle in the axisymmetric random flow dif- 
fers from that of a rod immersed in the flow resulting from the superposition of a uniform 



shear and of a short-correlated isotropic random component (|Puliafito fc Turitsvnll2005 



lTuritsvnl[2007h . In the presence of a strong mean shear, the tumbling dynamics of a rod 
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t 

Figure 5. A typical time evolution of $(t) in the tumbling regime (7 = 1, b/a = 4.75, 
c/a — 10, d/a — 0.9, T>' R /a — 10 -2 ); ■dit) has been computed by numerically integrating (|3.12[1 . 



consists of aperiodic transitions between two unstable states: the one aligned with the 
direction of the shear and that anti-aligned with it. When a fluctuation takes the rod 
away from the aligned or anti-aligned state and moves it into the unstable region of the 
flow the mean shear makes the rod flip. By contrast, in the axisymmetric case, large 
deviations of the orientation of the particle from the axis of symmetry of the flow do 
not necessarily result into sudden flips of the particle (figure [5]) . Simply, the orientation 
vector of the particle fluctuates randomly, but the orientations aligned and anti-aligned 
with the axis of symmetry of the flow are much more probable than the other orienta- 
tions. Thus, there can be excursions of $(t) from d « to d « ir/2 followed by a return 
to "d « (figure [5]). This behaviour would not be possible in the presence of a strong 
mean shear. 

The axisymmetric random flow and the laminar uniform shear differ in the dependence 
of the orientation dynamics on the geometrical shape of particles. In a uniform shear 
flow, the motion of a r od or of a dis c represents a degenerate case of the dynamics of 
non-spherical particles (Jeffcry 1922). Moreover, the dynamics is qua itativcly different 
for | 7 | < 1 and I7I > 1 ljBrethertonlll962l) . In the axisymmetric random flow, the dynamics 
of particles changes smoothly as a function of the shape coefficient. 

Finally, it is worth remarking that the function \ defined in (|4.5|) also determines the 
probability distribut ion of Jeffery's orbits in the presence of weak Brownian fluctuations 
( Leal fc Hinchlll97ll ). There does not seem to be, however, a simple relation between the 
orientation dynamics in the uniform shear and that in the axisymmetric random flow. 
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Appendix A. Rotary diffusion 

Rotary diffusion is introduced in the dynamics of N(t) by assuming that its effect is to 
produce isotropic Brownian fluctuations of the direction of N(t), but not of its length. 
Mathematically, this is obtained by adding a Laplacian term to the equation for f(n;t) 
which acts only on t he orientation of n. For a determini s tic flow, f(n; t) t hus satisfies (e.g. 
Leal fc Hinchlll97ll; iHinch fc Lea3ll972t lBrenner|[l973 iBird et al\\l9m iDoi fc EdwarcN 
19861) : 



3 t / = -V -[w(n, t)/] + ^V 2 /, 



(Al) 



where 



(A 2) 

with E(n) = I — nn/\n\' 2 and V„ = (8/3rii, d/dn2, 8/8713). The differential operator V 
is the angular part of the gradient or, equivalently, the restriction of the gradient to the 



V = r(n) • v„ 
12 and V n = (9/3711, d/dn,2, 8/8713) 



sphere of radius \n\. In (|A ljl 
The matrix IT satisfies: 



V = V • V and w(n, t) = n(t) ■ n — [n(t) : nn]n/ 



n x [V n • F(n)] = 0, 



r(n) 



(A3) 



By using properties (|A 3p and w(n,t) • n = 0, it is possible to rewrite (IA 1[) as follows: 



(A 4) 



at oui I oni ouj 

Equation (|2.3p is the Stratonovich stochastic differential equation associated with (IA4[) . 

It is worth noting that given that the orientation vector has unit length, /(«.; t) must 
take the form f(n;t) — F(n;t)5(\n\ — 1) with dF/dn — 0. It would therefore be more 
natural to consider the evolution equation for F(n;t) instead of t hat for f(n\t)\ this 
is indeed the usual approach in the literature jLeal fc Hinchl Il971[ IHinch fc Leall Il972 



Brennerlll974t IBird et allll977t IDoi fc EdwardJll986l ). In the present context, however, 



it is easier to formulate the problem in R 3 and to move to angular variables afterwards. 



Appendix B. Positive semi-definiteness of the covariance tensor 

The fourth-order tensor f is the covariance of a Gaussian second-order tensor, and 
must therefore be positive semi-definite, i.e. 

u t] r lopq u pq ^Q (Bl) 



for all second-order tensors U. Inequality (|B 1|) c an be reinterpre ted within the theory of 
positive semi-definite second-order tensors (e.g. iMoakher 20081) . Consider an invertible 



map I which assigns to each pair of indices (i, j), 1 ^ i,j ^ 3, a single index ranging 
from 1 to 9. By means of the map t, U can be regarded as a 9-dimensional vector; likewise, 
r can be regarded as a symmetric 9x9 second-order tensor, whose symmetry follows 
from Fij vq = r pqi j . Accordingly, (|B 1[) can be rewritten in the following form: 

l^l(i,j)Ap,q)^9 

Inequality (|B 2[) is the definition of positive semi-definiteness for second-order tensors. 
The theory of such tensors says that a necessary and sufficient condition for an Hermi- 
tian second-order tensor to be positive semi-definite is that all the principal minors of 
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the tensor are non- negative ( Gantmacher 19591 p. 307). When this condition is applied 



to rtu,j)tlp,q) it yields the following inequalities: 




and 



vji = 4a - b + d ^ 0, 



VJ3 = 15a 



if a ^ 
if a < 



vj 2 = 4a + 6 + c + 6d ^ 0, 
d)(c + 5d) + 2a(2c + 13d) ^ 0. 



(B3) 



(B4) 
(B5) 



In this paper, the above inequalities are assumed to hold strictly; (|2.15[) follows from this 
assumption. 

Some useful inequalities can be derived from (|2 . 1 5[) . Here, 7 ^ (the case of a spherical 
particle is indeed trivial), T>' R — (it is easily seen that a positive T>' R does not alter the 
inequalities below), and ^3 7^ (the case ^,3 — is treated separately in appendix |D|) . 

Firstly, A^{"&) and A^^tf) are positive for all values of a, b, c, d, and 7 satisfy- 
ing (|2.15p . Indeed, the quantities 2/xi and 2/ii + ^3 + ^4 are quadratic polynomials in 7, 
their discriminant is equal to —073 < 0, and for 7 = 1 they are equal to vj\ > and 
to Z02 > 0, respectively. Moreover, 2/ii + ^5 is a quadratic polynomial in 7 with dis- 
criminant equal to —(3a + d)(5a + d) < and takes the value 4a + d > for 7 = 1. 
Hence 



/ii>0, 2/xi + M5>0, 
for all a, 6, c, d, and 7, and consequently 



2/Wi + ^3 + ^4 > 



2/xx > 

lim j4 w (t?) = lim .4 W ($) = +i 



19— >7I 



and 



Am {■&) ^ m sin 2 (i9) + (2/ii + /x 3 ) sin 4 (i?) ^ (2/Xi + ^3 + Hi) sin 4 (i?) > 
> 2^1 + ^5 > 



(B6) 

(B7) 
(B8) 



(B9) 
(B10) 



for all < ■& < n. 

Secondly, ^ s t(i?) is bounded for all a, b, c, d, and 7 and for all ^ 1? ^ 7t. For Z\ > 0, 
this property is obvious. For A < 0, note that 



(>/=Z + A*4)(a 



M4j 



-M3-Pi(7)) 



(BH) 



and 



(\/^Z + M4 + 2^3)(\/^-^4-2/i3) =-M 3 P 2 (7), (B12) 
where -Pi (7) and ^2(7) are quadratic polynomials in 7 such that Pi(l) = iiui > 
and p2(l) = 4^72 > 0. Furthermore, the discriminants of Pi (7) and 1*2(7) are equal 
to — 16073 < 0. Therefore Pi (7) and Pa (7) are positive for all 7, and the products on the 
left-hand-sides of (|B lip and (|B 12j) have the same sign as —(13. 
Also observe that 



(\/^Z\ + ( U4 + 2 A i 3 )(\/^2\-/i4) =-M3 P 3 (7)-2 



(B13) 



where ^3(7) is a quadratic polynomial satisfying: Pf (7) + 4A = Pi (7)1*2 (7) > 0. 
Hence P 3 (7) - 2 v / ^3 ^ for all 7. Moreover, for 7 = 0, P 3 (0) - 2y/^A~ = P 3 (0) = 
vji + vj 2 + 2(a + d) > 0. As a result Pa(7) - 1\J-A > for all 7. The left-hand-side 
of (|B 13[) therefore has the same sign as —^3. 
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Three cases should now be distinguished: 

(a) ^3 < 0: in this case, \J — A ± ^ 4 > and \J— A — /14 — 2p% > (remember 
that A = 8/^1^3 — p\ with p,\ > 0), whence \J — A + /i 4 + 2^ 3 > (see (IB 12[) V Therefore 

2fi 3 sin 2 (#) ^ \/^Z\ 



A*4 



^4 + 2^3 > 0, 



(B 14a) 
(B146) 

~A - /i 4 < 



(B15a) 
(B156) 



^4 - 2/x 3 sin (1?) ^ V-^ - ^4 > 0. 

(6) ^3 > and /14 > 0: then \J~^A + /14 > and consequently 
(see (|B 11)) ). Hence 

\/ = 2i + fx 4 + 2(i 3 sm 2 (i9) ^ J~A + /Lt4 > 0, 
- /i 4 - 2/i 3 sin 2 (tf) s$ ^J~A - p A < 0. 

(c) /LA3 > and ^4 < 0: for these values of the parameters, (|B 12[) and ()B 13[) yield the 
following relations: 

^J~A + /i 4 + 2^ 3 sin 2 (i?) < ^J~A + ^4 + 2^ 3 < 0, (B 16a) 

\/~^A - /i 4 - 2[i 3 sm 2 (d) ^ ^J~A - p, 4 - 2p 3 > 0. (B 166) 

Inequalities (|B 14)) . ()B 151) . and (|B 16)) guarantee that if A < 0, the function is 
bounded for all a, 6, c, d, and 7. 

For Z\ = 8/ii/i3 — /i 2 = 0, /14 cannot be zero since \x\ > and /13 ^ 0. Moreover, ^3 
and hence c must be positive, in that p,\ > 0. Therefore, for ^4 > the function xW is 
bounded. The case /14 < requires a more detailed analysis. Z\ can be rewritten thus: 
A = c(6a + 5cf)7 2 (7 2 - p) with 

(26 + c + 5d) 2 -c(10a + c + 9d) 

P = ra l ■ ( B17 ) 

c(6a + 5a) 

(Note that, for 7 = 0, \i\ = (10a + c + 9d)/8 and hence 10a + c + 9d > for all a, c, d; 
this is consistent with the positivity of the variance of the components of the vorticity — 
see (|2.13[) .) Provided that p > 0, A vanishes for 7 = ±7* with 7* = yfp (the case 7 = 
is not considered here). Now note that both ^4 and 2p% + p,^ are quadratic polynomials 
in 7. If 7** = 1 + (26 + 5d)/c is positive, then 



|/i 4 >0 if 7 G (0,7**) 
\M4<0 if 7 ^[0,7**] 
If 7** < 0, then 

f/i 4 >0 if 7 e (-7**,0) 
\/i 4 <0 if 7 ^[-7**,0] 

Hence 

^4(2^3 

Let us now show that I7+I < |7+*| 



and 



and 



2M3 + M4<0 if 7 G (-7**,0) 
2p 3 + p 4 >0 if 7 i [-7**,0]. 

f2^ 3 + ^4<0 if 7 6(0,7**) 
}2^ 3 + /i 4 >0 if 7 ^[0,7**]. 



(-M4) > 
The quantity 7^ 

P 4 (c) 



V |7l < |7**l- 



7 2 is written: 



(B18) 



(B19) 



(B 20) 



(B21) 



(6a + 5d)c 2 ' 

where Pa{c) is a quadratic polynomial in c whose coefficients depend on a, 6, and d and 
such that P 4 (0) = (6a + 5d)(26 + 5d) 2 > and lim^oc P 4 (c)/c 2 = 4n?i > 0. Moreover, 
Pi(c) = if and only if c=c= -(26 + 5d)(12a - 26 + 5d)/(8roi), and for c = c: 

(26 + 5d) 2 [16nj 3 + (26 + 5d) 21 



Pi (?) 



16g7i 



> 0. 



(B22) 
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Therefore, Pi(c) > for all c > and 7^ > 7 2 for all c > 0. As a conclusion, if A = 0, 
then (|B 20|) holds, and for ^4 < 0: 

Hi + 2^3 sin 2 < ^4 + 2^3 < V < < 7t. (B 23) 

This result proves that xW is bounded also for A = and ^14 < 0. 



Appendix C. Stationary probability density function of the 
orientation angle 

Consider the functions: 



<P(d) = In 



Oo 



gg(g) 



dz 



and 







The stationary solution of (|3.7[) is written (see lRiskenlll989l p. 

3^(0) = ATe-'W - 5 S (i?)e-*W, 



(CI) 



(C2) 



(C3) 

where A/" and 5 are constants. As Awi^) and 2?,j(i?) are periodic, also e - *^- 1 is periodic. 
By contrast, g(-d) cannot be periodic given that A^$(^) > and hence g'("&) > 0. There- 
fore, #st(i?) satisfies (l3~T0)) if and only if S = 0. Equation (l4~2|) then follows from (|C 3[) 
with 5 = 0. 



Appendix D. The case ^3 = 

If /X3 = and 7^0 (i.e. c = 0), three cases should be distinguished. 
For /i 3 = 0, (1*4 7^ 0, and /Z2 ^ 3/^/2, the integral in (|4.2I) can be easily calculated by 
means of the transformation y = sin 2 (a;) to yield: 

VfetW =AA(2 Ml + M4 sin 2 7?)^~ i . (Dl) 



The stationary p.d.f. is bounded as a consequence of (|B 6j) . By examining the first deriva- 
tive of the above function, it can be shown that, depending on the value of the parameters, 
tpst has either three or five extrema in the interval 1? n. Therefore, the four regimes 
identified for ^ 7^ also describe the dynamics of the particle for fi 3 = 0. 
For /13 = /i4 = (i.e. c = 2b + 5d = 0), the stationary solution is: 



it(^) = A/"exp 



^sin 2 (tf) 
2/ii 



(D2) 



For 112/^1 > 0, ^st h as two minima in and 7t and one maximum in 7t/2, and hence the 
particle rotates in the plane orthogonal to A. For 1^2/^1 < 0, V'st has two maxima in 
and 7t and one minimum in 7t/2; therefore the particle tumbles between the direction 
parallel to A and that antiparallel to A. 

Finally, for fi 3 = and /12 = 3^/2 (i.e. c = 2b + (5 + "f)d = 0), (|D 1|) implies that the 
stationary statistics of orientations is isotropic: Vst($) = (47T)" 1 . 
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